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A new minimal model is constructed for the doped manganese oxides which exhibit colossal mag- 
netoresistance (CMR), involving a broad spin- majority conduction band as well as nearly localised 
spin-minority electron states. A simple mean field analysis yields a temperature-dependent hy- 
bridised band structure with suppressed carrier weight at the Fermi level. Spin stiffness is complex, 
indicating strong spin wave damping. Further investigations are needed to verify the relevance of 
the proposed model. 

PACS numbers: 75.47.Gk, 75.30.Mb, 75.10.Lp, 75.47.Lx 
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The unusual properties of doped manganese oxides ex- 
hibiting colossal magnetoresistance (CMR) |l| are not yet 
understood theoretically, and the problem of formulat- 
ing a suitable microscopic model remains open. These 
properties include the CMR phenomenon and metal- 
insulator transition, which are, in turn, intimately related 
to the temperature-induced variation of the electron den- 
sity of states ("pseudogap" [2]) or of the effective carrier 
number Q. This implies that, alongside double exchange 
ferromagnetism, the effects of electron-electron interac- 
tion play a key role and, quite possibly, are responsible 
for CMR itself. Indeed, the bandstructure calculations 
[3] suggest that the on-site Hubbard repulsion U is the 
largest energy in the problem. Its effects are typically 
considered within the one- or two-orbital model whereby 
the strongly-correlated behaviour is induced by an inter- 
action between the two spin-majority electronic e g states 
on-site, or between the two spin-components of a single 
e g bandjlj]. 

It should be noted that the experimental data[6j in- 
dicate the presence of spin-minority electrons near the 
Fermi level even in the low-temperature ferromagnetic 
state. This agrees with the bandstructure calculations 
@; 0] > suggesting that a narrow spin-minority band lies 
close to the Fermi energy. Since both the localised ti g 
and itinerant e g states originate from the same d-shell of a 
Mn ion, and therefore are characterised by approximately 
the same value of Hund's rule splitting Jh, it is clear (see 
Fig. Q] o) that these spin-minority electrons populate the 
spin-down ti g (localised) states [J]. This is further corrob- 
orated by the studies of Lai_ x Ce x Mn03 (with 1 + x con- 
duction electrons per formula unit) [H , which show both 
the spin-minority character of the carriers and a large 
overall increase in the resistivity (in comparison with the 
usual 1 — x-electron case, e.g., Lai_ x Ca x Mn03), consis- 
tent with the extra x electrons going into the spin-down 
t2 g states. Yet, while "two-fluid" models involving itin- 
erant and localised states were suggested earlier by both 

ll|, the appealing pos- 
of these states having 
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FIG. 1: (a) Crystal field splits the d-electron states of a Mn 
ion into t2 S (solid) and e g (dashed); these split further due to 
a lattice distortion. Spin degeneracy is broken due to Hund's 
rule, and the "spin-plus" ti g electrons form the local spin 3/2. 
Chemical potential (dotted line) lies within an e 9 band, with 
a "spin-minus" t^ g level nearby, (b) A simplified model of 
Eq-HJ), with the spin quantisation axis fixed by a local spin S. 
(c) Relevant orbitals in the limit Jh — > oo (see text); adding 
an itinerant (localised) electron increases (reduces) the total 
on-site spin of S + 1/2 by 1/2. 



antialigned spins has not been addressed theoretically. 

Furthermore, we note that electronic properties of a 
model where the orbital degree of freedom is taken into 
account are strongly coupled to lattice dynamics via the 
Jahn- Teller effect [l(, which results in additional split- 
ting of both e g and ti g levels. However, the phenomenon 
of CMR occurs in a broad class of bulk systems [three- 
dimensional (3D) perovskites and quasi two-dimensional 
(quasi-2D) bilayered] and thin films of varying chemical 
composition, and is presumably always due to the same 
physical mechanism. It is therefore worthwhile to con- 
sider a minimal model with fewer orbitals, which still 
captures some of the important intra-atomic physics of a 
Mn ion, before pursuing more complicated (and probably 
more material-specific) options 13 [ . 
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In the present Letter, we introduce such a simplified 
description and proceed with a simple mean-field analy- 
sis. While finer theoretical tools are required to gain a 
fuller picture, qualitatively our results for electron disper- 
sion and magnetic properties appear very encouraging. 

We consider a model involving a large spin 5 and 
two conduction-electron orbitals (broadened and nearly- 
localised) at each Mn site (cf. Fig. Q] b): 
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Here, the operators Ci a {di a ) annihilate an e g fog) con- 
duction electron of spin a =^, 4- (in the laboratory frame) 
at a site i of a square (or simple cubic) lattice. Localised 
spins Si originate from the remaining two ti g electrons, 
hence, in reality S = 1. They interact via superexchange 
J and are also coupled to the spins of conduction elec- 
trons on-site, = ^Yl a ,p ^p( c ia°i^ + ^L^w) (where 
&ot.p are the Pauli matrices) via a strong ferromagnetic 
Hund's rule exchange Jh] the external magnetic field 
H is applied along the z-axis. Owing to the fact that 
the (electron) co-ordinate operator is not diagonal in the 
band index, there is a hybridisation V between the ti g 

and e g states; E^ is the (bare) energy of the tig elec- 
trons. Direct hopping between ti g states on different sites 
is assumed to be negligible, while the nearest-neighbour 
hopping t between the e g states will be used as an energy 
unit, t = 1. 

We construct the spin-wave expansion for the model 
(HJ), kee ping terms up to first order in 1/5. In the spirit 



of Ref. [14[, it is expedient to introduce a new basis of 



electron states on each site according to 
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(expressions for df and dj, are obtained by substituting 
g a fa)- Operators 54 and fa (54 and fa) cor- 
respond, respectively, to the e g and ti g electrons with 
a spin parallel (antiparallel) to the total spin on-site, 
% = Si + <7i, whose vibrations are annihilated by a 
Holstein-Primakoff magnon operator /?,. The Hund's 
rule term takes form 
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Finally, the Holstein-Primakoff operators a% of the 
original spins Si are expressed as 

n 9f9i + ff h 9f9t + ff ft - fyl - ft fl p 
a ~ P 7= ! ! — P- 
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We next substitute these expressions into the Hamil- 
tonian ([IJ and take the limit of Jh 00 while keeping 
Ed = E^ + Jh constant. Now, if the chemical potential, 
denoted fj, + U — Jh/%, lies within the spin-up e g band, 
then the spin-down e g band is completely empty, and 
the spin-up ti g band completely filled. This is precisely 
the case of interest to us, containing the effects of the 
Coulomb repulsion between the spin-up e g and spin-down 
tig electrons in the presence of a localised spin 5+1/2 
(Fig. [1] c). Henceforth, we drop all the terms containing 
the operators (4 and ff, and suppress the spin index of 
remaining fermion operators. 

Assuming the ferromagnetic ground state, the Hamil- 
tonian takes the form % e + H.^ + T-L m with the electronic 
and (magnon-assisted) hybridisation terms, 

= E - m) g\g% + (E d m) E fci + 
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Here, N is the number of lattice sites, 



cos k T 



cos ky(— cos k z ) is the tight-binding dispersion law in two 
(three) dimensions, and Rj is the radius- vector of site j. 
It is assumed that the site basis is more appropriate for 
describing the narrow-band fermions fj, easily localised 
by fluctuations or disorder. The term "H m , which is of or- 
der 1/5, contains Zeeman electron energy shifts, magnon 
dynamics, and double-exchange band-narrowing effects: 
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Here, momentum-conserving summation is denoted by 
E', and d is the dimensionality of the system (2 or 3). 

The electronic term, H e , is the familiar Falikov- 
Kimball model. The rich physics contained therein (IBj 
crucially depends on the presence (and nature) of the 
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inter-band hybridisation. The form of our Tih, Eq. ([3]), 
is dictated by spin conservation: electron transfer be- 
tween the two opposite-spin bands must be accompa- 
nied by magnon creation or annihilation, pi or f3^. Such 

transfers require misalignment of spins Si on neighbour- 
ing sites, hence the hybridisation matrix element van- 
ishes at q — > 0, underlining the importance of short- 
wavelength processes. The latter feature appears promis- 
ing in the context of CMR compounds, where the unusual 
short-range correlations are reflected in the electronic and 
magnetic 3, 17 1 properties. Importantly, magnetic field 
H affects the carriers both via double exchange mecha- 
nism and by changing the energy difference between lo- 
calised and itinerant states [see Eq. Q]. 

We shall be interested in the regime characterised by 
non-zero values of both fermion occupancies n 9 and w 
(the latter assumed independent on R i ) 1 



In addition, there also arises an off-diagonal average, 
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(here and below, we omit the site indexes of operators fj 
once the j-dependent exponent has been factored out). 

In order to clarify the basic physics contained in our 
model, Eqs. (HHU, we will now proceed with a mean- 
field analysis of it. Here, we focus on the simplest self- 
consistent scheme, allowing only for average values ((SE]) 
and for a non-zero magnon occupancy. While actual va- 
lidity of this approach is probably restricted to the in- 
termediate temperature range (on the scale of the Curie 
temperature) and moderate values of U (see below) , it of- 
fers important guidance for future investigations. Mean 
field equations can be found in a standard way by de- 
coupling the equations of motion for the appropriate re- 
tarded Green's functions, expressing the latter as 
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Here, the magnon occupancy is 

(ni and are expressed in a similar way). Hartree 
energies of magnons and those of e g and ti g electrons 
read 

k 

h = h + Unf - t§ + ^ E^ e £f? - e £)> ( 13 ) 
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^- = ^+^E«S+2s. (14) 

The last terms in Eqs. (|12j|13[) contain the familiar double 
exchange physics (ferromagnetic contribution to the spin- 
wave spectrum and the magnon-induced conduction band 
narrowing, respectively) . 

In analysing Eqs. ([7 HTU)) we consider a 2D system; the 
3D case can be expected to be similar. We replace all 
factors (tk — e k-q) 2 with their average values over the 
isoenergetic surfaces e = and e = eq. We arrive 
at a system of mean field equations for /i, n^, and four 
quantities J$(e)(2 + e)vde,fN(e)tvde, fj\f(e)(v 2 ) £ vde, 
and JjV r (e)(4— e 2 )vde, where i/(e) and (i> 2 ) e are the tight 
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binding density of states and average velocity square at a 
given energy e. Solving these equations numerically, we 
observe that: 

(i) The localised band is broadened, and a temperature- 
dependent gap A (a new small energy scale) opens in 
the spectrum of itinerant electrons (fermions g£). Quasi- 
particle weight of itinerant electrons decreases when the 
energy approaches the gap from either side. The Fermi 
level lies below the gap, and the quasiparticle weight at 
the Fermi surface is strongly suppressed (Fig. [2] a and 
b). This behaviour, which is already reminiscent of a T- 
dependent (pseudo)-gap found experimentally @, will be 
further modified in a more exact treatment (e.g., a finite 
relaxation time will arise in the second order in 1/S). 
Eqs. ((TJ- HJ imply that in the present model, these spec- 
tral features are directly controlled by spin dynamics. 

(ii) With only the spin-majority electrons contributing to 
the spin stiffness D, which in 2D or 3D is given by 



dependence, especially in the region around and above 
Tc- In particular, this should lead to an unconventional 
longitudinal spin dynamics [20j and to a renormalisation 
of the Curie- Weiss constant (cf. Ref. 21]). In prin- 



ciple, the value of (T) should also be accessible more 
directly via muon spin rotation [22[ and NMR [23| mea- 
surements. We suggest that these methods (combined 
with electron spectroscopy) should be used to measure 
the value of (T). Its temperature dependence, especially 
if it correlates with (magneto)transport properties, would 
imply that a successful theoretical description of CMR 
compounds should indeed include spin-minority localised 
electrons. 

It is a pleasure to thank R. Berkovits, G. Goobes, A. 
Kanigel, K. A. Kikoin, E. M. Kogan, B. D. Laikhtman, 
D. Orgad, and S. Satpathy for discussions. This work 
was supported by the Israeli Absorption Ministry. 
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the latter is suppressed in comparison with the usual 
double exchange case (Fig. [2] c). An unusual feature 
of the present model is that the spin-flip continuum ex- 
tends down to zero energy and momentum [19|. As a 
result, D also develops an imaginary part. This im plie s 
strong magnon damping, as observed experimentally 17 1 . 
Magnon damping proportional to q 2 is not usually ex- 
pected in a ferromagnet 17ll24| and means that the spin- 
diffusion coefficient acquires a real (dissipative) part. 

As for the diffusive central peak found in the inelas- 
tic neutron scattering [l^ . we expect it to arise once the 
magnon-assisted diffusive motion of t2 g electrons (ne- 
glected here) is taken into account. Experimentally, 
strongly damped magnons, central peak, and pseudo- 
gap in the density of states @ (or optical Drude weight 
reduction^) are the key generic features of CMR man- 
ganates at the intermediate-to-high temperatures below 

A relatively small value of U used in Fig. [2] is due to 
the reduced stability region for mean field solutions with 
both n/ and n g different from zero. This reduction is 
an expected artifact of a simplistic mean field approach, 
mirroring, e.g., the greatly enhanced mean field stability 
of ferromagnetism in the Hubbard model. This situa- 
tion calls for further theoretical investigation, combining 
advanced mean-field schemes with numerical methods. 

These future treatments will also have to address the 
issue of ferro- to paramagnetic transition and a possibil- 
ity of charge ordering. We expect that any transition will 
be accompanied by a change of electron distribution be- 
tween the two bands, thus changing the magnitude of the 
net spin % on-site. Experimentally, the relevant quantity 
is the average total spin (T) = (3 + n 9 — n?)/2 of a Mn 
ion, which should show temperature and magnetic field 
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FIG. 2: Mean field results for a 2D system with Ed = —0.35, V = 0.2, U = 1.2, and J = 0. (a) Schematic view of the itinerant 
fermion dispersion (solid lines); dotted lines show the unhybridised e"k and Ed- (b): Temperature dependence of quasiparticle 
weight Z at the Fermi level (solid line) and the gap, A(T) (dashed line), for the electron density 1 — x = 0.7. (c) Doping 
dependence of real (solid line) and imaginary (dashed) parts of spin stiffness D at T = 0.0015. Dotted line shows spin stiffness 
for a usual double exchange model. t-2 g occupancy increases from n; — 0.25 at 1 — a; = 0.65 to nf = 0.43 at 1 — x = 0.73. 
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